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The Kondo effect in a Luttinger Uquid is studied using the renormaUzation group method. By 
renormalizing the boson fields, scaling equations to the second order for an arbitrary Luttinger 
interaction are obtained. For the ferromagnetic Kondo coupling, a spin bound state(triplet) can be 
realized without invoking a nearest neighbor spin interaction in agreement with the recent Bethe 
ansatz calculation. The scaling theory in the presence of the scalar potential shows that there is 
no interplay between the magnetic and non-magnetic interaction. Also a study on the crossover 
behavior of the Kondo temperature between the exponential and the power law type is presented. 



I. INTRODUCTION 



The Kondo problem in one-dimensional quantum system has attracted substantial interest in connection with the 
recent rapid development of the nanofabrication technology. In one-dimension, interacting electron system is described 
by the Tomonaga-Luttinger liquid theory [Q-^, whose low- energy excitations are not quasiparticles but collective 
charge and spin density fluctuations. The magnetic impurity efl'ect in this non-Fermi liquid was first studied by Lee 
and Toner using a scaling analysis on the kink gas action They obtained the scaling equation to the first order 
of coupling constant in the weak coupling regime. Furusaki and Nagaosa(FN) extended the study and obtained a set 
of scaling equation up to the second order using the poor man's scaling theory Q. In their study, FN proposed an 
interesting conjecture that even ferromagnetic Kondo impurity as well as antiferromagnetic one would be completely 
screened. In the antiferromagnetic coupling case, the coupling constant flows to a strong coupling regime and, thus, 
the magnetic impurity and a conduction electron form a singlet. But, in the ferromagnetic case, situation is not clear 
whether it flows into a strong coupling regime or not. In order to clarify the situation, FN considered coupling of 
impurity spins not only with the same site electrons but also with the nearest neighbor conduction electrons. In this 
picture, the impurity and three conduction electrons form a singlet composite. However, a recent Bethe ansatz(BA) 
result by the Wang and Voit showed that spins form a triplet for the ferromagnetic coupling case ||^ . These conflicting 
results call for a more detailed renormalization group(RG) analysis to clarify the physical situation. 

In this paper, we carry out an RG analysis which goes beyond the poor man's scaling scheme in the presence of the 
Kondo interaction. We bosonize the Kondo interaction term using the Abelian bosonization. Carrying out a full RG 
calculation for the second order cumulant, we obtain a set of general scaling equations valid for an arbitrary strength 
of the Luttinger interaction. We show that the present RG calculation conflrms the calculation drawn from the BA 
calculation using an open boundary condition. 

Another recent BA study of the same problem by Wang and Eckern showed that there is a competition between 
the Kondo coupling and the impurity potential When the impurity potential is dominant, the system is shown to 
flow to a weak coupling fixed point. However, when the magnetic interaction is dominant, a spin complex is shown 
to be formed. In order to investigate this problem, we bosonize the potential scattering term and perform a scaling 
calculation. The result shows that the magnetic interaction and the potential scattering do not interplay in the scaling 
procedure and fiow independently. 

In the first order scaling, the Kondo temperature in the Luttinger can be calculated analytically to show a power-law 
behavior in contrast to the exponential one of the conventional Kondo model Q. However, it is not possible to show 
the crossover from the Fermi to the Luttinger liquid in the first order. In this paper, we study the crossover behavior 
of the Kondo temperature from the exponential to the power-law type as a function of the Luttinger interaction, by 
solving the RG equations. 



II. BOSONIZATION OF THE KONDO HAMILTONIAN AND THE PARTITION FUNCTION 



Introducing the two phase fields, 

in ^ g-a|p|/2-»p2. 
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the one-dimensional model Hamiltonian with the forward electron-electron scattering is written by 



dx 

where v = p or a, ) means right(left) going mode, and the parameter Vi, and 7],^ are given by 
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Here, 52 and 54 are the dimensionless coupling parameters The fields (/>y and tt^ satisfy the canonical boson 
commutation relation; 

[0i,(x),7r,,/(x')] ^ iS^^'S{x - x'), (4) 

[0i/(a;),6',,/(x')] = i-(5yi./sign(a;' - a;), d^{x) ^ it '!r^{y)dy. 

^ J -00 

Alternatively, this Hamiltonian can also be obtained directly from the fermion representation of the one dimensional 
interacting electron system using 

ir{kp — j-\x 

where 77^3 is the Majorana fermion operator which satisfies the following relations |^; 

[rirsTrjr's']='iSrr'5ss' (6) 

If a single magnetic impurity is introduced in one-dimensional interacting electron system, the Kondo interaction 
term is given by 

Hk = JS ■ s(0) 

= J.5,s,(0) + i Ji(5+s_(0) + 5-s+(0)), (7) 

where s= ^J2rr'aa' i^la-'^crcr'''Pr'a' and s± = Sa; ± isy. Using the relations of the Majorana fermions, Eq.(6), and the 
bosonization formula of the fermion operators, we obtain the bosonized Kondo Hamiltonian [ p^ , 

Hk^ JS-s 

= -^{JzFad^V2(j)^[Q) + 2iJ,BV+]V-] sin(%/2<?!>p(0)) cos(V20<,(O))) 

/TTQ! 

+ ^e-^'''''''\j^FV+}V+l cos{V2M0)) - U^BV-^V+l sin(\/20p(O))). 
The partition function of the system at temperature T = 1//3 is 

Z = j D^pD(l)„D6pD6,e~^, (9) 
5*= j '^^ j d-'^i^o + Lk), 



Lo= E {idr4>u{x,T)'n:^{x,T) + ^{vv'q^TT'^Til 

v=p,a \ \ 



Vv ( d^uix.T) 



LA'=i?A'(</'.(O,T),0,(O,T)), 



2N 



r\i, \ dx 

where the integration is over the bosonic fields ip^ix, t) and Oy{x, r) with imaginary time r running from to /3. 
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III. RENORMALIZATION ANALYSIS 



First, we divide the phase fields into slow and fast mode: 

|CJ„|</J 



Lo e 



(10) 



{u)e 



M<l^nl<A 

The average over the fast mode of the partition function is carried out, using the cumulant expansion. 



(11) 



where Zq = f D(j)^D9^e~^° , f indicates an average over the fast mode which will be omitted hereafter, and < • • • >o 
represents an average over Zq. 

First, we consider the first order forward longitudinal scattering term, 



(12) 



The second term in the right side vanishes because it is an average of an odd function. We, thus, conclude that 
SJzF = in the first order. The backward longitudinal scattering part is scaled as follows, 



dTj^BV+Tn-^ Sm{V2(f)ps{T) C0s{V2(f>crsiT)). 



x) 

In the above, we utilized the fact that the charge and spin degrees are separated. The rescaling procedure, 

^p+Vf7 -I 



Jzsip) = 



Jzb{\) 



gives 



where fj, = X + SX and 51 



Jz 
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(5 In A. Thus, we have in the first order 

dJzF 



dl 

dJzB 

dl 
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J. 
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The scaling equations for the other scattering terms can be similarly obtained; 

f )] ^±F, 

|)1 Ae. 



dJ±F 


1 - 




dl 




dJ±B 


1 - 




dl 





(17) 



These equations are in agreement with those of Lee and Toner 

The second order cumulant is given by — ^(< S]^ > — < Sk >^), where < Sk >^ term is to eliminate unconnected 
diagrams. We consider one of the JzbJ±f terms which is given by 
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J dr J dr' -^^J;,B'2iil+]r]^]Jj_FV+lV+i: 

(< sin(y20p(r))cos(\/20,(r))e^^^'("''cos(y20,(r')) > (18) 
- < sin(V20p(T))cos(y20,(T)) >< e^*^-("')cos(V2(/),(T')) >) . 

In order to evaluate this expression, we need the two point correlation function |TT- 14 

G(x, r) = < (f>{x, r)(/.(0, 0) > (19) 



27r J 2tt -^uj^ + ii(72 

G{t)= G(0,t) 

^Koifir) for At » 1 



^Inf forAT«l, 

where Kq is the modified Bessel function of the second kind. G(t) decays exponentially with the renormalized lattice 
spacing i, and decreases logarithmically for small At. Thus, we regard G(t) short ranged and, thus, expand the 

cosine terms around zero. It is known that the higher harmonics cos {2V24>{Q)) is irrelevant and {{d<l>{T)/dT)\r=Q)'^ 
terms which is the most relevant term in the expansion can also be shown irrelevant by the power counting 
Short time cut off, tq ^ a/vp, merely introduces an overall constant which does not affect the flow of the parameter. 
Therefore Eq.(18) is reduced to 

/ (2^5*--'-^-(-'«''-'' iijl"" - ^^"'"^■»(^*,(^))- (») 

Collecting other terms of the second order cumulant and rescaling as before, we have for the transverse backward 
part, 

dT^ir,+^r^^^ri,J,BJ±F^e^'''^^^^ sin(V20p,(T)), (21) 
Zna A 

which renormalize the coupling constant J±b of the Kondo term. Similarly, we have for the transverse forward 
scattering part which renormalizes the coupling constant J±f, 

dT^r^+^r^+^ripJ,BJ±B^e^'''<-^^ cos(V2^4t)). (22) 

The same scaling process on the S- terms gives the same renormalization for both the J±b and J±f- 

However, the scaling process involving the descendant field terms, JzFO:dxV24'a{T), is somewhat different. One 
example of such a term is the second order cumulant for the transverse forward scattering, which is given by 

- < adxV2(t>„{T) >< e'V2«<'(^') cos{V2(t>„{T')) >) (23) 
- ^'^+rJ.FJ±FV+lV+1 < a5,V2</)„/(T)e'^«'(^')cos(V2</)„(T')) > ■ 



(27ra)2 



Here, we note that 



< aa,y2</),/(T)e'^«'(^')e^^'^"(^') > 

= lini 1 A < „g..y2*./(x,r)g^y2e„(O,r')g»y20.(O,r') > | ^ (24) 

Using e^^^ — e^e^ e^^^^'^\ the above expression takes a form, 

-{2dxG^,,{x,T)~2dx<<f>,f{x,T)e,f{Q,Q) >)|,^oe-^'''/^°^°)-^^*/^°^"^ (25) 
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In order to calculate dx < 4>a-f{x,T)9c,f {0,0) >, we use the relations between density field, (f>a-, and current field, 9^ 



i d(t>a 


_ 06^ 


Vcr d{v^T) 


dx 








dx 



(26) 



Then, we have 



adx < 4>afix,T)9^f{0,0) > = at]^ 



d 



d{v„T) 
IdX 
2T' 



<9,f{x,T)9,f (0,0) > 



(27) 



using Eq.(19) with jj ^ rj and the short-ranged nature of the correlation function G(r). Substituting this result into 
Eq.(23), we obtain 



1 

27ra 
1 

2'Ka 



(^J.FJ^By) e^^''"(^)sin(V20p(r)), 



(28) 



for the renormalization of J±f and J±b- 

The longitudinal scattering parameter is scaled by consecutive transverse scattering processes. The forward scat- 
tering part which contains J^p is 



S+S- 
(27ra): 



;JlpV+lV+^V+mi \< e^^'-'^)cos(V2^.(r))e-^^«"(^')cos(V20.(r')) 



> 



< 



> 



The typical relevant term is given by 



gV2ie„(r)g-V2ie„(r')g-V2i0„(T)g%/2i0„(T') 



(29) 



(30) 



Other terms containing higher harmonics, e^^'^"^'^) or e2V2i0ff(r) ^^.^^ irrelevant pT| , p^ . Separating the fields into the 
fast and slow mode, we average on the fast mode to obtain a short range correlation and expand the slow mode to 
obtain 



1 + iV2dr9crs{T)\T=oT + i\/29^ (/f^rs (t) | ,-=0^ - 



(31) 



We note that dr9a-s{T) gives the longitudinal scattering contribution, dx4>as, through Eq.(26), while dr4>iys term is 
canceled by other terms. Including the contributions from the S-S-^- term, we obtain 



dT 



1 I ( I , 7/, 



2TTa 2ttvf T]<j \2?7ct 2 

2TTa 2'KVp rjcr \2r]cr 2 



J 



dX 
A 



adxV2(t>crsiT) 



(32) 



for the J^g contribution to the J^f renormalization. Consecutive transverse scattering also renormalizes the longi- 
tudinal backward term, JzB, similarly; 



dr;^2i?7+t77_t- Ty^r J±_f J±b^ sin(V20p(r)) cos(V2(^o-(t)). 

Zna liTVp A 



(33) 



The last term to be considered in the second order cumulant expansion is the S^Jfp term. Scaling of this term, 
however, does not contribute to the magnetic interaction and, also, does not produce any relevant term for non- 
magnetic interaction. Collecting the results together, we now have the scaling equations to the second order; 
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dJzF 1/1/1 \ ,2 , W 1 \ j2 

^ = ^ U ''^^ + ^ ) 

dJzB /. Vrr -Vr,\ ^ 1 



dJ±F 
dl ' 

dJ±B 




2t],J^fJ±b, (34) 
{JzfJ^f + 'r]pJzBJ±B), 
{JzfJ±B + V<rJzBJ±F)- 



dl 

These equations can be simplified assuming the SU(2) symmetry for the conduction electrons {rj^ = 1 

Jl = {pqJ±f) + [-^^ j iPoJ±B) 

dipoJzB) _(l-Vp, ^^^j^^^ ^ 2{poJ^f){p,J±b) (35) 

{P^Jzf){PoJ^f) + 'npiPoJzB){PoJ±B) 
{PqJ±b) + {PqJzf){PqJ^b) + iPoJzB){pQJ±F), 



dl 

d{poJ±F) _ 

dl 

d{poJ±B) _ f ^-Vp 

dl 



where po{— \/2'KVp) is proportional to the density of state. As we see later, for ferromagnetic J, the triplet ground 
state can be formed due to the interacton parameter -qp in the second order term. In ID correlated electron system, 
the electron-electron interaction induces a short range magnetic ordering which, in turn, introduces a molecular field 
on the magnetic impurity [ p^ . Therefore, the impurity and a electron form a triplet aligned to the molecular field, 
thus, resulting the broken local SU(2) symmetry. These equations are naturally reduced to the scaling equations of 
the conventional Kondo model when there is no Luttinger interaction i.e. t]^ = 1 and r]p = 1. 

So far, no assumption has been made on the strength of 77p(= [(1 — g/7rwi?)/(l + g/T^VF)]^)- Here, g is the strength 
of scattering between the right-going and the left-going mode and the scattering within the same mode is neglected. 
In the small g regime, the above scaling equations are reduced as follows, 

d{poJzF) _ T ^2 , /„ T \2 ^„ T ^2 



(PoJ^f) + (poJ^b) ~ (poJ^b) 

Zttvf 

(PoJzb) + 2{poJ±F)ipoJ±B) (36) 



dl 2t:vf 
djpoJzB) _ 9 
dl 2'hvf 



d{poJ±F) ^ {pqJzf){PoJ±f) + {PoJzB)ipoJ±B) - ipoJzB){PoJ±B) 

dl Zttvf 

d{poJ±B) 



dl 2'kvf 



(poJj-b) + {poJzf){poJj-b) + {poJzb){poJ±f)- 



We note that the above result is in agreement with the FN's poor man's scaling result except for the symmetry 
breaking terms involving the g parameter. In fact, the extra terms in our RG equations correspond to the next higher 
order terms which were neglected in the FN approach. 



IV. THE GROUND STATE 



The scaling equations yield two strong coupling fixed points; {Jf,Jb) = (oo,oo) and (cx3,— oo). The first fixed 
point governs the antiferromagnetic regime, which gives the singlet as the ground state. The second one corresponds 
to the ferromagnetic coupling. However it is not clear whether this fixed point corresponds to a singlet or to a triplet 
state. We have calculated flows of the coupling constant for several values of the Luttinger interaction parameter, 
?7, for the ferromagnetic fixed point. For finite rj other than unity, Jb/Jf < —1 and Jb + Jf flows to — oo, whereas 
Jb + Jf is equal to zero when rj = l(Fig. 1). When the coupling constant J grows and becomes large, the Kondo 
coupling term becomes dominant in Eq.(35). Therefore, the Luttinger interaction becomes irrelevant and we can 
treat the impurity spin as a classical spin with the magnitude S = 1/2. 

The asymtotic Hamiltonian which represents this situation can be written as 
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H ^^ekclf^^Crka + ^ Jrrv4/caC^'fcV', (37) 

rka rk,r'k',(T 

where Jno- = ^220- = ctJf and Ji2a = ^210- = ctJb- The Luttinger interaction effect is included in the coupling 
constant through the renormalization process. The corresponding Green's functions are given as follows ; 

Gn(e) = G22(e) = Go(e) + Go(e)- , f'^^ G,{el 

Gi2(6) = G2i(e) = Go(6) ^ j"'^'' Go{e). (38) 

1 - j(7{Jf + JB)Go(e) 

For T] — 1{Jf + Jb — 0), there exist no bound state since the Green's functions do not have poles. This case has the 
same properties as the three-dimensional ferromagnetic Kondo coupling case. However, if the Luttinger interaction is 
turned on(?7 < 1), Jp + Jb flows to —00. The Green's functions have a pole for a = +1 suggesting that a conduction 
electron becomes bounded to form a triplet state in the ground state for the ferromagnetic Kondo exchange. 



V. THE SCALAR POTENTIAL SCATTERING EFFECT 

In a more realistic model, a magnetic impurity generates scattering due to an elastic potential, uj | |To[ |. Recently, 
Wang and his coworkers have carried out Bethe ansatz calculations on the one-dimensional Kondo problem . They 
showed that for an attractive potential scattering, the bound state always appears, whereas for a repulsive potential 
scattering, there exists a bound state for \J\ > 4w, but no bound state for | J| < Auj. 

The scaling theory for the scalar potential in the first order is similar to the magnetic interaction. The bosonized 
term is 

-^{2LUFad,V2<j)p{0) + iiujBV+lV-1 cos(\/20p(O)) sin(\/2</),(0))), (39) 
Zira 

and we have the scaling equations 

duJF _ „ diOB _ I - rip 



dl ' dl 



LOB, (40) 



showing that the backscattering contribution is relevant. In the second order, the candidates for renormalizing the 
forward potential scattering are contributions from the terms like ujp or w^. But, we find that tOp yields only a 
constant by a simple calculation, 

j^^AuW JdT J dr' < d.,V2(bp, (t)9,\/20p, (t') >, (41) 

which clearly has no slow mode. Similarly, the ujg term generates terms such as cos(2\/2(/>o-(''')), and cos(2\/2(/>p(t)), 
which are irrelevant as discussed above. Also, we can show that J^p, J^g, JzfJzB, wfwb and Jz±,fbujfb are also 
irrelevant through the same procedure as in ujp and ujg. Thus, the potential scattering is not scaled at the second 
order cumulant and there is no interplay between the magnetic and the non-magnetic interaction. Actually, this 
situation is not unexpected, since we know that the nearest neighbor coulomb interaction in the Hubbard model does 
not affect the magnetic properties jist . 

In the case of an impurity spin, S=l/2, the energy of a single electron coupled with the impurity is such that the 
ferromagnetic coupling energy is |J and the scalar potential is uj. It means that if 4uj > \ J\ initially, the potential 
energy grows infinitely along the scaling process so that no bound state can be formed. If Auj < \ J\ initially, a spin 
triplet state is formed. Under the open boundary condition, this corresponds to a bound state of a spin | complex in 
agreement with the Bethe ansatz results Here, we note that, in the present treatment, it is not necessary to invoke 
the nearest neighbor spin interaction as done by FN. We believe that the short-ranged nature of the spin interaction 
makes the FN scenario unlikely, although it cannot be ruled out completely. 
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VI. THE KONDO TEMPERATURE 



In the conventional three-dimensional Kondo model, the Kondo temperature is given by Tk = De 2j ^ which 
originates from the scaling equation dJ/dlnD = — 2J^, where D is the bandwidth. In their previous study, Lee and 

2 

Toner showed that the scaling relation, dJB/dlnD — —^-^^Jb, gives a power law Kondo temperature, Tk — DJg""^ 
However, the crossover behavior from the Fermi liquid to the Luttinger one has not been studied. In order to address 
this question, we consider the scaling to second order for the forward and the backward scattering simultaneously. 

From the RG flow, we get the coupling constants as functions of D, i.e. Jf{D) and Jb{D), where Jp and Jb are 
the scaled couplings and D is the scaled bandwidth cutoff. The Kondo temperature is an invariant energy scale in the 
scaling procedure and, thus, can be expressed as Tk = Df{J), where f{J) becomes exponential or power law type 
depending on the limiting case. For some initial value of J, Tk remains constant through the scaling procedure. In 
such a case, /(J) becomes proportional to '^Jyj^ thus yielding the Kondo temperature. 

Tk for the backward scattering part is given in Fig. 2. It can be clearly seen that Tk is exponential for a weak 
Luttinger interaction and becomes power law type as the Luttinger interaction increases. From inset, we observe that 
the linear slope of linear plot(r/ = 1) is -1/2 in agreement with Tk = De~2j. 



VII. SUMMARY 



In summary, we have studied the Kondo effect in a Luttinger Liquid in presence of a scalar potential. We have 
obtained the scaling equations up to the second order for an arbitrary Luttinger interaction strength by renormalizing 
the boson fields. 

The ferromagnetic fixed point is studied using an asymptotic Hamiltonian. It is shown that a triplet bound state 
can be formed in agreement with the recent Bethe ansatz calculation without invoking the nearest neighbor spin 
interaction. The Luttinger interaction induce the triplet state to break the local SU(2) symmetry about the impurity 
spin contrary to the result obtained by Furusaki and Nagaosa. The magnetic interaction and the potential scattering 
do not interplay and the triplet state is sustained for a weak scalar potential, \J\ > The Kondo temperature 
for arbitrary strength of the Luttinger interaction is calculated. The result shows a clear crossover behavior from an 
exponential to a power law type. 
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Figures 

Fig. 1 The coupling constant Jf,b for several Luttinger interaction strength r]. The magnitude of slope and the 
value oi\JF + Jb\ become larger as the interaction strength grows larger. The slope of the dotted line is -1. 

Fig. 2 The Kondo temperature as a function of the backward scattering coupling constant, Js, for several Luttinger 
interaction strength 77, where is the Kondo temperature for each r] aX Jb / D = Q.l. The linear to nonlinear crossover 
is clearly shown. Note the log- log scale of the graph. The inset shows the Kondo temperature as a function of J^^. 
The linear curve for 77 = 1 clearly shows the Tk = De~^ behavior. 
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